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Orientational Melting in Carbon Nanotube Ropes 
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Using Monte Carlo simulations, we investigate the possibility of an orientational melting transition 
within a "rope" of (10, 10) carbon nanotubes. When twisting nanotubes bundle up during the 
synthesis, orientational dislocations or twistons arise from the competition between the anisotropic 
inter-tube interactions, which tend to align neighboring tubes, and the torsion rigidity that tends to 
keep individual tubes straight. We map the energetics of a rope containing twistons onto a lattice 
gas model and find that the onset of a free "diffusion" of twistons, corresponding to orientational 
melting, occurs at Toii/^IBO K. 
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Since their first successful synthesis in bulk quantity 
PI , "ropes" of single- wall carbon nanotubes have been in 
the spotlight of nanotube research. Recent experimental 
data indicate that carbon nanotube ropes exhibit an un- 
usual temperature dependence of conductivity |0] , mag- 
netoresistance ||^, and thermoelectric power [Q. Several 
physical phenomena have been suggested to cause the 
intriguing temperature dependence of conductivity be- 
havior, such as twistons g, orientational melting [|6|,[7|, 
weak localization |g[ , and Kondo effect g] . The opening 
of a pseudogap near Ep in a bundle composed of (10,10) 
tubes has been postulated to result from breaking the 
Dioh tube symmetry by the triangular lattice |^,p . 

Nanotubes |o) and Cgo "buckyball" molecules [|l]jl| 
are similar inasfar as their interaction is weakly attrac- 
tive and nearly isotropic when condensing to a solid. 
The small anisotropy of the Cgo intermolecular potential 
drives the solid to an orientationally ordered simple-cubic 
lattice with four molecules per unit cell at low tempera- 
tures ||l|l. Only at T>249 K does the Ceo solid undergo a 
transition to a face-centered cubic lattice. As confirmed 
by ^^C nuclear magnetic resonance |14 , this is an order- 
disorder phase transition, with Cgo molecules spinning 
freely and thus becoming equivalent above 249 K. Un- 
like the Ceo solid, very little is known about the equilib- 
rium structure of bundled nanotubes beyond the fact that 
they form a triangular lattice ||l|. In particular, nothing 
is known about the equilibrium orientation of the tubes 
within a rope. More intriguing still is the possibility of 
an orientational melting transition associated with the 
onset of orientational disorder within a rope. 

In the following, we calculate the potential energy sur- 
face and the orientational order of straight and twisted 
tubes within a rope. We further postulate that realis- 
tic nanotube ropes contain orientational dislocations that 
have been frozen in as the tubes formed ropes at a finite 
temperature, the same way as dislocations are known to 
form in crystals. We map the energetics of a rope with 
dislocations onto a lattice gas model and find that the 



onset of orientational disorder, corresponding to a free 
axial diffusion of twistons, should occur at Tom^160 K. 

In order to determine the orientational order and the 
rotational motion of tubes in a rope, we need to describe 
both the inter-tube interaction and the torsional strain 
within the individual tubes. Since the anisotropic part 
of the inter-tube interaction in a rope is weak and lo- 
cal, it can be well described by pairwise inter-tube {not 
inter- atomic) interactions. To describe the interaction 
between coaxial tubes as a function of their orientation, 
we use the parametrized linear combination of atomic 
orbitals (LCAO) formalism with parameters determined 
by ab initio results for simpler structures [|l5[ . This tech- 
nique has been successfully used to explain superconduc- 
tivity arising from inter-ball interactions in the doped 
Ceo solid jl^] , the opening of a pseudogap near Ep in a 
(10,10) nanotube rope |^,0, and the opening and clos- 
ing of four pseudogaps during the librational motion of 
a (5,5)@(10,10) double-wall tube (l|]. The total energy 
functional consists of a nonlocal band structure energy 
term and of pairwise interatomic interactions describ- 
ing both the closed-shell repulsion and the long-range 
van der Waals attraction. A smooth cutoff function |l^] 
has been implemented to keep the total energy continu- 
ous as the neighbor topology changes while tubes rotate. 
This energy functional correctly reproduces the exfolia- 
tion energy, the interlayer distance and the C33 modulus 
of hexagonal (AB) graphite, as well as the energy barrier 
for interlayer sliding, corresponding to the energy differ- 
ence between AB and AA stacked graphite. 

To determine the interaction between a pair of aligned 
(10,10) nanotubes, we first define the orientational an- 
gles ifi and ip2 for these tubes by the azimuthal angle of 
the center of a particular bond with respect to the con- 
nection line between adjacent nanotube axes, as shown 
in Fig. |y(a). For each ((/?!, (^2) pair, we calculate the 
inter-tube interaction using a fine mesh of 800 k-points 
sampling the one-dimensional irreducible Brillouin zone. 
Due to the high symmetry of the system, the interaction 
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FIG. 1. (a) Definition of the orientational angles (fi and 
ifi2 of two aligned nanotubes. (b) Interaction energy per atom 
between two (10,10) carbon nanotubes at equilibrium distance 
as a function of the tube orientations, (c) Top view of the 
equilibrium structure of bundled interacting nanotubes, with 
a two-tube unit cell, (d) Torsion energy per atom within an 
individual (10,10) nanotube. 



energy AE[ipi,ip2) is periodic in ipi and ip2, with a pe- 
riod A(p = 36°. With the simple variable transformation 
9i = (fi + ip2 and O2 = ifi — (p2, the interaction energy 
AE can be simply expanded in harmonic functions of 9i 
and 62, with the same period A9 — 36°. The resulting 
potential energy surface, displayed in Fig. |l](b), shows a 
maximum corrugation of only AEmax~0.5 meV/atom. 

The equilibrium geometry of bundled (10,10) nan- 
otubes can be found by optimizing the total energy E 
with respect to the orientations of all individual tubes 
p9| . Due to the high level of orientational frustration 
in a triangular lattice of (10,10) nanotubes with Dio/i 
symmetry, the potential energy surface E(ipi,ip2, ■ ■ ■) is 
very complex. We determine the global minimum of 
E{ipi,(p2, ■ ■ ■) by applying the Metropolis Monte Carlo 
algorithm to an infinite system of straight nanotubes in 
a (2x2), (4x4), ... , (8x8) superlattice, with unit cells 
containing between 4 and 64 nanotubes. Independent of 
the unit cell size, we find that the global energy min- 
imum corresponds to a two-dimensional oblique lattice 
with two tubes per unit cell, shown in Fig. ^(c). In equi- 
librium, the orientations ipi, ip2 of the two tubes within 
this unit cell satisfy the condition 9i = (pi+(p2 = 12° and 
O2 = ^i-<y32 ~ ±9.71° within the range -18°<e' < +18°. 

Even though the energy barrier AE'^O.S meV per atom 
for a free rotation in a pair of tubes, shown in Fig. 0(b), 
appears small, the barrier to rotate an entire tube, that is 
completely straight and rigid, is unsurmountable. In the 
following, we postulate that tube rotations in a rope are 
still possible in view of the finite, albeit large value of the 
tube torsion constant. To determine the torsional strain 
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FIG. 2. Temperature dependence of (a) the total energy 
IS.E and (b) the specific heat C per degree of freedom in model 
(mxm) superlattices of interacting (10,10) nanotubes with 
zero torsional rigidity. Each nanotube is modeled by a chain 
of torsionally decoupled, rigid segments ("disks") containing 
twenty atoms, each representing one orientational degree of 
freedom. In this model system, orientational melting occurs 
at T>roM~55 K. 



within an isolated, twisted (10, 10) tube, we combined 
the LCAO method mentioned above with the recursion 
technique ||2^] . This approach has been used successfully 
to describe the dynamics of fuUerene melting pi|, the 
growth of multiwall nanotubes 123] , and the dynamics of 
a "bucky-shuttle" memory device Q. Our calculations, 
shown in Fig. 0(d), suggest that the torsional energy is 
harmonic up to a strain of Aip / ALk,1° / K. Within this 
harmonic regime, the torsional energy per atom can be 
well represented by the expression AE = k(A</7/AL)^, 
with Kfti2.58xl0~^ meV- rad^/A^. Since the number of 
atoms in the tube is proportional to the total tube length 
ALtot, the total torsional energy of the tube is inversely 
proportional to the tube length for a given total twist 
angle Aiptot and hence becomes vanishingly small for a 
long tube. 

To study the possibility of orientational melting in a 
nanotube rope, we first consider an unrealistic model 
system of bundled (10, 10) tubes consisting of torsion- 
ally decoupled axial segments of 20 atoms, resembling 
rigid "rings" or "disks". The interaction between two 
adjacent disks in neighboring tubes is the value AE of 
Fig. pFb), multiplied by the number of atoms in the 
disk. Absence of axial coupling makes this system equiv- 
alent to a two-dimensional triangular lattice of disks with 
one (orientational) degree of freedom per disk. Results 
of Monte Carlo simulations of orientational melting in 
(2x2), (4x4), and (6x6) superlattices with 4, 16, and 
36 such tubes per unit cell, respectively, are shown in 
Fig. g. To ensure proper phase space sampling even at 
low temperatures, each data point represents an ensemble 
average taken over ^10^ Monte Carlo steps per degree of 
freedom. Results for the temperature dependence of the 
total energy per degree of freedom, shown in Fig. 0(a), 
suggest that an orientational melting transition should 
occur at Tom~55 K- This transition becomes more pro- 
nounced with increasing unit cell size in the superlat- 
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FIG. 3. Temperature dependence of (a) the total energy 
AE and (b) the specific heat C per degree of freedom in a 
realistic rope consisting of interacting tubes of finite rigid- 
ity, (c) Illustration of the microscopic process of orientational 
melting in an exaggerated perspective. The two views depict 
the position of orientational dislocations at difi^erent points in 
time at r = 500 K. The axial motion of the twistons, indi- 
cating orientational melting, is highlighted by the changing 
grey-shaded sections in four of the tubes. 



tice. The sharp peak in the temperature dependence of 
the corresponding specific heat data, shown in Fig. y(b), 
suggests this phase transition to be of first order. At low 
and at high temperatures, the specific heat approaches 
the classical value C = 0.5/cb. 

To address orientational melting in a realistic (10, 10) 
nanotube rope, we add the proper torsional coupling be- 
tween the disks, according with our results in Fig. 0(d). 
Monte Carlo simulations analogous to those described 
above showed no indication that initially straight nan- 
otubes would start to perform a "ratcheting twisting mo- 
tion" , corresponding to orientational melting in the rope, 
up to 4, 000 K when individual nanotubes should disin- 
tegrate structurally. 

The key to the understanding of orientational melt- 
ing is to consider the energetics and dynamics of orien- 
tational dislocations in the system. Consider two nan- 
otubes which, under synthesis conditions, show a total 
twist A(/Jto(>36° over the entire tube length, at the neg- 
ligible cost in total energy of wQ.l eV for a 100 ^m long 
tube. As these tubes bundle up during annealing, at 
least two orientationally aligned domains form within 
the tube pair, separated by an orientational dislocation. 



Similar twists of up to l°/nm, associated with such ori- 
entational dislocations, have been recently observed by 
electron difi^raction in nanotube ropes [g^ . Such solitons 
may move rather freely along the tube axes, but cannot 
be annihilated if the paired tubes are infinitely long, or 
if they bundle up to a double-torus [^ . 

The equilibrium geometry of these frozen-in twistons 
is given by the competition between the anisotropic 
inter-tube interactions, which tend to align neighboring 
tubes, and the torsion rigidity that tends to keep indi- 
vidual tubes straight. Their energetics and dynamics 
is well represented by mapping the Frenkel-Kontorova 
model onto the orientational degrees of freedom of a 
triangular lattice of chains. In bundled (10, 10) nan- 
otubes, the orientation (p{z) of individual tube layers 
near the dislocation is well reproduced by the function 
(p{z) — ipo + Aip (1 -|- exp ({z — zq)/w)) , with the typi- 
cal values A(^f«36° and w«30 A. Such dislocations, which 
can be either left- or right-handed, show a total twist of 
36° that extends over w250 A in the axial direction. The 
energy cost of Ai?w5 eV to create such a twiston, given 
by the energy difference between ropes containing one or 
no dislocation, is relatively high. This explains why no 
twistons or left- /right-handed twiston pairs were created 
in perfectly straight tubes even at high temperatures. 
Since the twistons cannot be created easily, the dynam- 
ics of the entire system is limited to a subspace of the 
configurational space, where the number of twistons on 
each tube is fixed. We notice that the axial motion of 
a twiston corresponds to a finite tube rotation in that 
given segment, and that the energetics of a rope con- 
taining twistons can be mapped onto a lattice gas model 
of twistons moving along individual tubes. The small 
potential energy barriers associated with an "up"- and 
a "down" -moving twiston passing each other in adjacent 
tubes depend on the orientation of the other surrounding 
tubes, and are evaluated by the total energy expression 
above. We correlate the onset of orientational melting in 
the rope with depinning and a completely free "diffusion" 
of twistons within the tubes. 

Results a Monte Carlo simulation for the orientational 
melting transition in a (10, 10) nanotube rope containing 
twistons are presented in Fig. ||. We consider a (6x6) 
superlattice of nanotubes with a fixed number of orien- 
tational dislocations on each tube. The energy of this 
system is given by the positions of these twistons within 
each periodically repeated unit cell containing 36 tubes 
with 5,000 layer segments discretizing the axial direction. 
This system with nominally 180,000 degrees of freedom 
is mapped onto a lattice gas of twistons in the following 
way. Out of the 36 tubes per unit cell, we select twelve 
non-adjacent tubes, each containing a single twiston that 
can move axially. The remaining tubes in the unit cell 
have eight such orientational dislocations frozen in. Their 
positions are equally spaced over the tube axes, but axi- 
ally offset in adjacent tubes. The factual impossibility to 



change the number of twistons, discussed above, is mim- 
icked by assuming the same handedness for ah twistons. 

The temperature dependence of the total energy of this 
classical system, shown in Fig. ||(a), shows a slope that 
is initially small close to T = 0, then becomes large 
within the temperature range K< T<500 K, and fi- 
nally becomes small. The corresponding specific heat 
data, shown in Fig. 0(b), begin with the classical value 
C = 0.5fcs at low temperatures and peak at Tqm^ISO K, 
the orientational melting transition. The nature of this 
transition is illustrated in Fig. |3|(c). The two snap-shots 
of the geometry, taken at T = 500 K, indicate that above 
T~ToM, frozen-in twistons become depinned and diffuse 
relatively freely along the tube axes, as indicated by the 
changing gray-shaded sections in the tubes. Since these 
twistons are important scattering centers for electrons, 
their depinning may significantly affect the transport in 
these one-dimensional systems at Tk^Tom that may play 
the role of T* in Ref. |]. 

In summary, using Monte Carlo simulations, we have 
investigated the possibility of an orientational melting 
transition within a "rope" of (10, 10) carbon nanotubes. 
We postulate that during the synthesis, as twisting nan- 
otubes bundle up, orientational dislocations or twistons 
arise from the competition between the anisotropic inter- 
tube interactions, which tend to align neighboring tubes, 
and the torsion rigidity that tends to keep individual 
tubes straight. We have mapped the energetics of a 
rope containing such twistons onto a lattice gas model 
and find that the onset of a free "diffusion" of twistons, 
corresponding to orientational melting, should occur at 
ToAf >160 K. 
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No. N00014-99-1-0252. 
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